Crystal structure and lattice dynamics of A1B 2 under pressure 
and implications for MgB 2 
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The effect of high pressures to 40 GPa on the crystal structure and lattice dynamics of AIB2 
was studied by synchrotron x-ray powder diffraction, Raman spectroscopy, and first-principles cal- 
culations. There are no indications for a pressure-induced structural phase transition. The Raman 
spectra of the metallic sample exhibit a well-defined peak near 980 cm -1 at GPa which can be at- 
tributed to the Raman-active E2 S zone-center phonon. Al deficiency of ~11% in AIB2, as indicated 
by the x-ray data, changes qualitatively the electronic structure, and there are indications that it 
may have a sizable effect on the pressure dependence of the E2 g phonon frequency. Similar changes 
of the pressure dependence of phonon frequencies, caused by non-stoichiometry, are proposed as an 
explanation for the unusually large variation of the pressure dependence of T c for different samples 
of MgB 2 . 

PACS numbers: PACS: 63.20.-e, 74.25. Kc, 78.30.-j, 71.15.Nc, 62.50.+p 



I. INTRODUCTION 

The AIB2 structure type and derivatives thereof are 
among the most frequently occurring-ones for intermetal- 
lic binary and ternary compounds.EH3 Transition metal 
diborides, belonging to this family, have been studied 
in some detail because of their potential application in 
electronic devicesa to overcome current problems of elec- 
tromigration, corrosion, and diffusion into the semicon- 
ductor substrate. The largest interest, however, has un- 
doubtedly received the recently discovered superconduc- 
tor MgB2 which also crystallizes in the simple AIB2 struc- 
ture depicted in Fig. |l|. 

AlB2-type compounds have not been studied systemat- 
ically at high pressures. There seem to exist no confirmed 
reports on pressure-induced structural phase transitions 
in metal diborides. There are, however, structural studies 
of the rare-earth metal digaUid£aXidGa2, HoGa2, ErGa2, 
and TmGa2 at high pressures. QoErQ With the lighter rare 
earth elements La — Er the digallides crystallize in the 
AIB2 structure at ambient pressure. In essence, all of the 
rare-earth metal digallides studied so far show a transi- 
tion to the UHg2 structure which is isotypic to AIB2 but 
with a lower c/a ratio. 

The superconductor MgB 2 was studied at high pres- 
sures with jCftg^d-.ta its superconducting, transition. 
temperatureMiMBBM .-crystal structure, HeNO'EI 
and lattice dynamics. t3'L3 The pressure dependence 
of T c could well be explained in the fra mework of 
phonon-mediated, i.e. BCS, superconductivity EJoEj An 
isostructural transition near 30 GPa was reportedE3 but. 
could not be reproduced in another study up to 40 GPa.Ea 
Much higher pressures may be necessary to induce tran- 
sitions possibly towards the UHg 2 structure. 

We study here the effect of hydrostatic pressure on 
the crystal structure and lattice dynamics of AIB2. Syn- 




FIG. f: Crystal structure of AIB2. The boron atoms form 
honeycomb layers. Al atoms are located at the centers of 
hexagonal prisms that are formed by the B sheets. 



chrotron x-ray powder diffraction and Raman scattering 
experiments are complemented by first-principles calcu- 
lations. The present high-pressure study was in part 
motivated by the question whether some of the un- 
usual physical propect^ of MgB2 - such as the siz- 



able anharmonicityllj'E 



the remarkably large calcu- 
lated pressure dependence^ of the E2 g phonon - are spe- 
cific to MgB 2 or whether they are characteristic of other 
AlB2-type metal diborides as well. It aims at providing 
high-pressure structural and lattice dynamical informa- 
tion for comparison with corresponding data on MgB 2 . 
Finally, Al deficienc#,a»pears to be hardly avoidable in 
the growth of A1B2.E3E3 We show that it has significant 
effect on the electronic structure of AIB2 and may also 
influence its lattice dynamics. We will discuss possible 
consequences of metal deficiency for the superconductor 
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MgB2 where this issue is also of relevance. 



II. EXPERIMENTS 
A. Experimental Details 

The structural properties of AIB2 under pressure were 
studied up to 40 GPa by monochromatic (A = 0.3738 A) 
x-ray powder diffraction at the European Synchrotron 
Radiation Facility (ESRF Grenoble, beamline ID30). 
Commercially available AIB2 powder (Alfa Aesar, 99%) 
was placed in a diamond anvil cell (DAC) for pres- 
sure generation. Nitrogen was employed as a pres- 
sure medium to provide nearly hydrostatic conditions. 
Diffraction patterns were recorded on image plates and 
then integrated^!! to yield intensity vs. 26 diagrams. 

Raman spectra of A1B 2 up to 25 GPa (DAC, 
4:f mcthanol/cthanol mixture as a pressure medium) 
were excited at 633 nm utilizing a long-distance micro- 
scope objective. They were recorded in back-scattering 
geometry using a single-grating spectrometer with a 
multi-channel CCD detector and a holographic notch fil- 
ter for suppression of the laser line (Dilor Labram) . For 
the Raman experiments the DAC was equipped with syn- 
thetic diamonds (Sumitomo type 11a) which emit only 
minimal luminescence. In all experiments pressures were 
measured with the ruby luminescence method.Ej 



B. X-ray diffraction under pressure 

Figure |^ shows x-ray diffraction patterns of AIB2 for in- 
creasing pressures up to 40 GPa. The diagrams evidence 
small amounts of Al metal as a secondary phase. At pres- 
sures above 2 GPa additional reflections are observed due 
to various phases of solid nitrogen. There are no indica- 
tions for a pressure-induced structural phase transition 
in A1B 2 up to 40 GPa. 

Lattice parameters as a function of pressure were deter- 
mined from Rietveld-type fits of the diffraction diagrams. 
The compressibility of A1B 2 is moderately anisotropic as 
illustrated in Fig. ||(a) with the softer direction being 
parallel to the c axis. Up to 40 GPa, compression along 
c is 47% larger than along a. The c/a ratio decreases 
from 1.083 (0 GPa) to 1.060 at 40 GPa. From the lat- 
tice parameters we determine the unit cell volume as a 
function of pressure as shown in Fig. ||(b). The data are 
well represented by the Murnaghan relatione^ V(P) = 
V [(B'/B )P+ I]- 1 / 8 '. With Va = 25.4734(5) A 3 fixed 
at the value determined from the zero-pressure data we 
obtain by least-squares fitting the bulk modulus Bq and 
its pressure derivative B' at zero pressure as listed in 
Table | 

Rietveld refinements of the crystal structure with the 
Al site occupation as a free parameter indicate that the 
sample studied here has an Al deficiency of ~11%. This 
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FIG. 2: X-ray powder diffraction diagrams of AIB2 at various 
pressures (T = 298 K). Al' marks a peak due to fcc-Al. Addi- 
tional peaks appear at pressures above 2 GPa due to various 
phases of solid nitrogen. 
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FIG. 3: (a) Experimental lattice parameters of AIB2 as a func- 
tion of pressure, normalized to their respective zero-pressure 
values, (b) Experimental (symbols and solid line) and calcu- 
lated (dashed line) pressure-volume relations. The lines are 
given by the Murnaghan equation of state with the parame- 
ters listed in Table |. 



is illustrated in Fig. ^ where difference curves for refine- 
ments of AI1.00B2 and AI0.89B2 are shown together with 
the experimental diffraction pattern. A Stephens peak 
profileEj was used and a common isotropic thermal pa- 
rameter for Al and B was optimized. The weighted profile 
R value (without background) reduces from R wp = 8.5% 
for AI1.00B2 to R wp = 6.9% for A1 . 8 9B 2 
of substantial Al defici 
density measuremeatsE 
diffraction results. E3 



This indication 
' in A1B 2 is in agreement with 
and recent single-crystal x-ray 
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TABLE I: Structural parameters of AIB2 and MgB2: zero-pressure volume Vq and lattice constants ao,co; bulk modulus Bo, 
and its pressure derivative B' at zero pressure. Variation of the c/a ratio with pressure is described by the quadratic polynomial 
c/a = co/ao + aP + (3P 2 with the coefficients a, j3 as listed below. 





Vo 

(A 3 ) 
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(A) 


CO 
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(GPa) 
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a P 
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A1B 2 , Exp. (300 K) 
A1B 2 , Calc. (DFT/C 


25.473(1) 
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FIG. 4: X-ray powder diffraction diagram of AIB2 at GPa 
and 293 K and difference curves (A = J cxp — J ca i c ) for refine- 
ments with and without Al-deficiency. Markers show peak 
positions due to AIB2 and fcc-Al. 
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FIG. 5: Raman spectra of AIB2 at ambient conditions 
recorded on two crystallites and at a sample spot where no 
crystallites were discernible ('powder'). 



C. Raman spectra at ambient conditions 

AIB3-. has one Raman-active zone-center phonon 
mode£3 It is an in-plane vibration of the B atoms with 
E2 g symmetry, where neighboring B atoms move out of 
phaseJij The powder sample we investigated contained 
some small shiny crystallites up to ~10 (im in size. Fig- 
ure [| shows Raman spectra recorded on two different 
crystallites and at a sample spot where no crystallites 
were discernible with an optical microscope. Besides a 
Lorentzian-shaped peak (FWHM of 40-50 ^em^ 1 ) which 
was attributed to the Ei g mode previously^ we observe 
an additional step-like feature at the lower-energy side of 
the main peak. It is clearly visible in the single-crystal 
spectra and reduces to a weak shoulder in the powder 
spectrum. The pcpsider spectrum resembles that reported 
by Bohncn et al.S3 The peak position of the Ei g mode 
of the two crystallites differs by 15 cm -1 (uo — 973 and 
988 cm -1 ). In the powder spectrum the main peak oc- 
curs at an even lower energy of 952 cm -1 . The step-like 



feature in the spectra of the crystallites shifts by the same 
amount as the main peak indicating that it is intrinsic 
to AIB2. It appears likely that it is related to a peak 
in the calculated phonon density of statesa which exists 
slightly below the energy of the E^g mode. 



D. Raman spectra under pressure 

Raman spectra of AIB2 were recorded for increasing 
pressures up to 25 GPa. Since the sample is somewhat 
deformed when pressure is applied and there seems to 
be some inhomogeneity of the powder, several spectra 
were collected at different locations of the sample. The 
spectra shown in Fig. ^(a) result from averaging about 
five spectra recorded at different spots that were selected 
for a narrow Ei g peak and low background. The zero- 
pressure frequency of the Ei g mode in these averaged di- 
agrams amounts to 981(1) cm -1 . Two additional peaks 
near 875 and 1025 cm -1 (at GPa) are due to the 
methanol/ethanol pressure medium. 
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FIG. 6: (a) Raman spectra of AIB2 recorded for increasing 
pressures of 0-25 GPa (T = 300 K). Peaks marked by as- 
terisks are due to the methanol/ethanol pressure medium, 
(b) Experimental and calculated energies of the E^ a phonon 
mode. Experimental data are represented by circles and a 
solid line, results of the GGA/frozen-phonon and LDA/linear- 
response-theory calculations by a dash-dotted and dashed 
line, respectively. The lines represent quadratic functions fit- 
ted to the data. 



are carried out with a plane- wave cutoff energy of 23.6 
Ry, and the Brillouin zone sampling is based on a T- 
centered 16x16x16 uniform mesh, yielding 270 to 1170 
fc-points in the irreducible wedge of the Brillouin zone, 
depending on the symmetry of the lattice (equilibrium 
structure or crystal with the displacements of the E2 S (r) 
phonon). As the system is metallic, the fc-space integra- 
tion with the incompletely filled orbitals uses the tetra- 
hedron methodO with BlochPs corrections .cj 

Phonon frequencies were independently verified us- 
ing the linear respopse theoryOE-3 as implemented in 
the ABINIT packageeS within the local density approx- 
imation (LDA)c3tD to the DFT. We used Hartwigsen- 
Goedecker-Hutter pseudopotentialscS treating Al(3s 2 3p) 
and B(2s 2 2p) levels as valence states, a plane- wave cut- 
off of 60 Ry, and a 12 x 12 x 12 fc-point mesh. A Gaussian 
smearing with a broadening parameter of 0.04 Ry was ap- 
plied to improve fc-point sampling in the special points 
method. The structural parameters were optimized for 
each volume/pressure such that the stress tensor compo- 
nents <j xx and <j zz agreed within 10~ 3 GPa. 



The E2 g phonon frequency increases continuously with 
increasing pressure, but with slightly decreasing slope 
dtu/dP. The zero-pressure mode-Griineisen parameter 
amounts to 70 = 1.22(3) (based on B Q = 170(1) GPa). 
Up to 10 GPa the peak width (FWHM) is essentially con- 
stant at 44 ±2 cm -1 while at higher pressures it increases 
to ~90 cm -1 at 25 GPa. This is most likely related to the 
solidification of the pressure medium near 10 GPa. From 
the Raman data there are no indications for a structural 
phase transition. 



III. CALCULATIONS 
A. Theoretical Method 

The theoretical methods employed here were described, 
in detail in the context of our recent work on MgB2.c2l 
In summary, first-principles calculations of the total en- 
ergy E t ot of the solid are the basis for the determination 
of the equation of state (EOS) and of the phonon fre- 
quencies. The energy is computed within the density 
functional theory (DFT) using a plane-wave basis and 
pseudopotentials . 

For part of the actual calculations (equation of state 
and phonon frequencies in the, fmrmrphonon approach) 
we employed the VASP codesEjO'L£lE3 within the gen- 
eralized gradient approximation (GGA)Ell The ultra- 
soft Vanderbilt-typej)seudopotentialsEa were supplied by 
Kresse and HafnerO The pseudopotential for Al treats 
explicitly three valence electrons (3s 2 3p) ; no semi-core 
states are included. A 'harder' variant of the potential 
was chosen and the nonlinear core correctionEJ is applied 
in order to improve the transferability. The calculations 



B. Structural Properties 

For fourteen unit cell volumes in the range 21.8-26.2 
A 3 we calculated (with the VASP codes) the total en- 
ergy E{V) for different c/a values, thus determining the 
optimized structures [a(U), c{V)] and minimized total 
energies. The latter data were fitted by the Murnaghan 
relation for ^(U)^ 



E(V) = E 



BqVq 
B' 



V 

To" 



(V/Vo) 1 



-B* 



B' 



B' - 1 



(1) 



which provided the static equilibrium volume Vq as well 
as the bulk modulus Bq and its pressure derivative B' 
at zero pressure. The structural parameters are summa- 
rized in Table [jj. The pressure-volume relation, shown 
in Fig. ^|(b), compares well with the experimental data. 
The calculated equilibrium volume is 0.4% larger than 
the volume measured at 300 K. The individual lattice 
parameters are —0.3% off the experiment for an, +0.9% 
for co, and consequently +1.2% for co/ao- The variation 
of c/a under pressure can be represented by a quadratic 
function and its coefficients are given in Table ||. 



C. Phonon frequencies 

Energies of the Ei g phonon mode as a function of pres- 
sure were initially calculated in the frozen-phonon (FP) 
approach using the VASP codes and the GGA. For each 
phonon and pressure, the atoms are given six different 
displacements ranging from u/a = 0.005 to 0.035 and the 
calculated energy values E(u) are fitted with a quartic 
polynomial. The resulting harmonic phonon frequencies 
and their respective pressure dependencies are listed in 
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Table [IT] together with the experimental data for thp-^g 
mode and our previous FP-GGA results for MgB 2 B As 
illustrated in Fig. ^(b) the calculated zero-pressure fre- 
quency is only ~1% lower than the experimental value, 
which is in the typical range for the difference between 
experiment and calculations of this type. Regarding the 
pressure dependence of the phonon frequency, however, 
we find an unusually large deviation of the theoretical 
results from the experiment. From the variation of the 
calculated phonon frequency with volume we obtain a 
zero-pressure mode-Gruneisen parameter of 70 = 1.00 
whereas the experimental value amounts to 70 = 1.22(3). 

Since the deviation of the theoretical results from the 
experimental data is larger than usual we performed 
a second calculation of the E 2g frequency under pres- 
sure using a rather different approach, namely linear- 
response-theory (LRT) in the LDA with the ABINIT 
package. Consistent with the common overestimation of 
bond strengths in the LDA we obtain here a somewhat 
larger zero-pressure phonon frequency, such that the two 
theoretical values bracket the experimental data at zero 
pressure. Over the whole pressure range the LDA/LRT 
calculation gives phonon frequencies which are consis- 
tently ~2.4% larger than the corresponding GGA/FP 
results. Consequently, we obtain essentially the same 
pressure dependence dlnw/dP for the two calculations. 
The deviation of the calculations from the experimental 
data is far beyond the typical uncertainty of such com- 
putations. The important difference between theory and 
experiment could be that the former is based on the ideal 
stoichiometry AI1B2 whereas the real sample is Al defi- 
cient. 

For comparison with previous calculations for MgB2 we 
have also calculated the frequency to of the B\ g phonon 
in AIB2 (out-plane motion of the boron atoms) Jl3 see 
Table ||. The pressure dependence is characterized by 
a mode-Griineisen parameter 70 = 1.06. In case of a 
constant, i.e., pressure-independent j(V) = 70 the re- 
lation u(V) = UoiV/Vo)" 10 holds. Thus, for both the 
E 2g and the B\ g mode in AIB2 with 7 w 1 there is a 
nearly inverse-proportional relation between the phonon 
frequency and volume. This is quite different from the 
situation in MgB2 where we do not only have the very 
large 70 = 2.5 for the E 2g mode as noted before but also 
a rather small 70 = 0.6 for the B\ g phonon (Table |J). 

The frozen phonon calculations also yield information 
on the anharmonicity of the phonon modes-as described 
in the context of our MgB2 calculations J2_3 In essence, 
the variation of the total energy with atomic displace- 
ment, Etot(u), can be represented by a polynomial where 
the ratio of the quartic to squared quadratic coefficients 
04/02 is a measure of anharmonicity. In the harmonic 
limit 04 = 0. For the E 2g mode in AIB2 we obtain 
|a4/o2| < 0.01 eV _1 which is about three orders of mag- 
nitude lower than the corresponding values for MgB2 of 
a^/a\ = 4-8 eV _1 (see Refs jl9| , p3| ). Small anharmonici- 
ties are calculated for the B\ g modes of both compounds: 
a A /a 2 2 = 0.27 for A1B 2 and o 4 /ol = -0.05 for MgB 2 . 



IV. DISCUSSION 

A. Structural stability 

In our x-ray diffraction and Raman experiments we do 
not find any indication for a structural phase transition 
or modulation of tha structure. Group-subgroup sym- 
metry considerations^ indicate numerous possible distor- 
tions of the aristotype AIB2 most of which are realized in 
intermetallic compounds at ambient pressure. Pressure- 
induced structural phase transitions of XY2 intermetal- 
lic compounds have not been studied systematically. In 
the context of the-AlB2 structure a number, of rare-earth 
metal digallidesJffllj'Ll KHg2,c3 and LaCu2E3 were investi- 
gated at high pressures. From the available data a transi- 
tion between the structure types AIB2 and UHg2 appears 
as a typical route. The KHg 2 (CeCu2) structure type 
may occur as an intermediate phase. AIB2 and UHg2 are 
isopointal structures, distinguished only by their c/a ra- 
tios. Two clearly separated groups of compounds of AIB2 
and UHg2 type are observed when plotting the c/a ra- 
tio versus ratiapij the metallic radii of XY2 intermetallic 
compounds.BOH The AIB2 and UHg 2 type branches are 
characterized by c/a ratios of 0.95-1.20 and 0.60-0.85, 
respectively. The compound A1B 2 with c/a — 1.083 (at 
GPa) is located near the center of the former branch. 
In the pressure range to 40 GPa explored here it de- 
creases only to 1.060. Pressures well above 1 Mbar may 
therefore be needed for a possible transition towards the 
UHg 2 structure. At lower pressures a transition involv- 
ing a buckling of the boron honeycomb layers may oc- 
cur which CQuJd lead to phases of the CeCu 2 , CeCd2, or 
Caln 2 typeiffl 

B. Raman spectra of AIB2 vs. MgB2 

Raman spectroscopy is commonly applied to semicon- 
ductors and insulators but only to a much smaller ex- 
tend to metals. It is essentially the group of elemen- 
tal hep metals that has been studied systematically, al- 
ready in the late 1960s at ambient pressureEil and more 
recently at high pressures (see e.g. Refs.|54|,|55] and refer- 
ences therein). It may therefore be attributed to a lack 
of reference data that in case of MgB2 the observation 
of a very broad Raman feature (FWHM of ~300 cm -1 ) 
near 600 cm -1 lead to a still unresolved controversy ovgij 
the origin of this peak. It has initially been attributedE2l 
to the Raman-active E 2g mode which immediately raises 
the question of the large linewidth. The large peak width 
has been related to both strong electron phonon cou- 
pling and to structural disorder. The latter now ap- 
pears less likely because Raman spectra of MgB 2 pow- 
ders are quite [Similar to those of recently available small 
single crystalsEEttJ which are presumably less disordered. 
High-pressure Raman experimentsEj have cast doubt on 
the assignment to the E 2g phonon. They revealed a 
double-peak structure with peaks at 603(6) cm -1 and 
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TABLE II: Pressure and volume dependences of selected phonon frequencies of AIB2 and MgB2. The zero pressure frequency 
coo and the linear and quadratic pressure coefficients were obtained by least square fits of u)(P) = uio + ct ■ P + f3 ■ P 2 to the data. 
The theoretical mode Griineisen parameters 70 (at the theoretical equilibrium volume) are derived from a similar quadratic 
expression for u(V). P has been obtained from V through the calculated P—V relation (see text). The experimental mode 
Griineisen parameter is determined from lo{P) and the experimental P—V relation. 



Compound 


Mode 


Wo 


a 





70 






(cm- 1 ) 


(cm^/CPa) 


(cm-VGPa 2 ) 




A1B 2 


E 2 g 


969 


5.425 


-0.0235 


1.00 




E 2g (exp.) 


981(1) 


7.027 


-0.0365 


1.22(3) 




Big 


490 


2.968 


-0.0087 


1.06 


MgB 2 


E 2 g 


535 


8.974 


-0.0780 


2.5 




Big 


695 


3.065 


-0.0190 


0.6 



750(20) cm -1 . Neither of the two peaks could be at- 
tributed to the E 2 g mode because of severe deviations 
from calculated phonon frequencies in terms of zero- 
pressure frequencies and/or the pressure dependences. 

The present Raman data of AIB2 show that it is pos- 
sible to obtain Raman spectra with a well-defined E 2g 
peak from metallic samples of the AIB2 structure. The 
difficulties encountered in case of MgB2 are therefore not 
likely related to the metallicity of the sample nor intrinsic 
to the structure type. Crystallinity also appears to have 
only a small effect on the Raman spectrum as the Ra- 
man linewidth of the AIB2 powder sample is comparable 
to that of the AIB2 crystallites. 

There are two properties of MgB2 with regard to 
phonons which make it distinct from AIB2. First, the 
whole E2 g phonon branch along the T-A direction in the 
Brillouin zone ex hibits v ery strong electron-phonon cou- 
pling in MgB2.cj'E§E3E3 Second, the E2 g (T) mode shows 
pronounced anharmonicityoO Both electron-phonon 
and phonon-phonon interaction decrease the phopDn life- 
time and hence increase the phonon linewidth.EJ They 
are therefore the most likely causes for the absence of a 
well-defined E 2g Raman peak in MgB 2 . 
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FIG. 7: Calculated electronic bandstructure and density of 
states (DOS) of AIB2 at ambient pressure. Energies are given 
with respect to the Fermi energy Ef- The DOS is given in 
units of (states/eV/formula unit). 



C. Metal deficiency in A1B 2 

Our x-ray diffraction data indicate an Al deficiency of 
11% in AIB2 in accord witb^-g¥evious density measure- 
ments and chemical analysisEaEias well as recent single- 
crystal x-ray diffraction results.t3 Although the change 
of the Ei g phonon frequency at different sample spots 
suggests that there is some variation of the Al content, 
there is no indication that growth of aluminum diboride 
in the composition AI1.0B2 is possible. The occurrence 
of substantial metal deficiency appears to be common to 
many (transition) metal diboridesJ£3 

In case of AIB2 this metal deficiency has important in- 
fluence on the electronic structure. A comparison of the 
calculated bandstructures of AIB2 and MgB2 shows that 
the relative ordering and dispersions of the bands near 
the Fermi level are very similar. The difference between 
AIB2 and MgB2 can largely be treated in a rigid-band 



picture with a higher band filling for AIB2 . It is therefore 
justified to discuss the observed Al deficiency of ~10% 
in a rigid-band picture, too. .Figure [?] shows the calcu- 
lated electronic bandstructureS and density of states of 
A1B 2 . In case of the stoichiometric compound the 7r-type 
bands derived from boron p z states and the bands with 
Al-s character (near T) are partially filled and give rise 
to the metallic state. 10% deficiency of Al removes 0.3 
valence electrons and consequently lowers Ep by 0.9 eV 
which leads to a complete depletion of the band with Al-s 
character. In other words, the experimentally observed 
Al deficiency is expected to lead to a qualitative change 
in the electronic structure compared to the ideal case of 
AI1B2. In terms of the Fermi surface this change means 
a removal of the electron pocket around T. 

The measured mode-Griineisen parameter of the 
AI0.9B2 sample is ~20% larger than calculated for 
AI1.0B2, a deviation which is far beyond the typical 
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uncertainty of such computations. The zero-pressure 
phonon frequency, on the other hand, seems to be hardly 
affected. The apparent difference between experiment 
and calculation is the Al deficiency of the sample which 
was not taken into account in the theory. We tenta- 
tively attribute the discrepancy between the experimen- 
tal and calculated mode-Griineisen parameters to the 
Al-deficiency-induced electronic changes discussed above. 
However, it cannot be excluded at this point that other 
effects are also at work. A more detailed analysis of the 
effect of Al deficiency on the electronic and lattice dy- 
namical properties of AIB2 is, however, beyond the scope 
of this work. 



D. Metal deficiency in MgB2 

Mg deficiency in MgB 2 was often indicated by the oc- 
currence of MgO as a secondary phase in MgB2 samples 
that were grown from a molar 1:2 mixture of Mg and B. 
Variations of T c and even more of the pressure depen- 
dence of T c for different samples appeared to be related 
to non-stoichiomctry of the material. The correlation 
between composition, structural parameters and T c has 
been established in an experiment by fndenbom et aL.Lj 
By diffusion of Mg into a boron cylinder they produced a 
sample with a composition changing gradually between 
Mgi.oB2 and Mgo.sB2. With decreasing Mg content the 
lattice parameter c increased by 0.003 A (0.1%) and T c 
increased fromJ17.2 to 39.0 K. 

Tissen et alHi furthermore pointed out a correlation 
between the zero-pressure critical temperature T c .o and 
the pressure derivative dT c /dP: As T Ci0 increases from 
37.3 to 39.2 K for various samples, dT c idP changes from 
—2.0 to —1.1 K/GPa. Monteverde et ala discussed a sim- 
ilar observation on a smaller number of samples in terms 
of the electronic band structure and band-filling effects 
related to the Mg non-stoichiometry. On the other hand, 
the pressure-induced changes of the electronic density of 
states calculated for MgB2 are too small to account for 
the observed decrease of T c under compression E3 It is 
rather the increase in the relevant phonon frequencies 
which provides the main contribution to the pressure de- 
pendence of T c £-jE2lEJ It would therefore be rather sur- 
prising if electronic density effects - i.e., the electronic 
density N{Ep) at the Fermi level - were responsible for 
the large sensitivity of dT c / dP on Mg non-stoichiometry. 

The indications that Al deficiency in AIB2 may affect 
the pressure dependence of the Ei g phonon frequency 
hints at an alternative possible explanation for the large 
sensitivity of dT c /dP on Mg deficiencjz-|in MgB2. It 
was first pointed out by Yilderim et alr3 that the E2 g 
phonon mode in MgB2 exhibits a very large anharmonic- 
ity. Boeri et alE3 showed theoretically that this effect 
arises in MgB2 because here the Fermi level Ep is located 
only ~0.5 eV below the top of a bands of the equilibrium 
structure. The lattice distortion of the E% g mode induces 
a splitting of these a bands large enough that the lower 



split-off band sinks completely below Ep£!3£!B This does 
not happen in AIB2 and graphite, and anharmonicity is 
indeed negligible. 

It is also noteworthy that the Ei g mode in AIB2 is 
much higher in energy than the B\ g phonon whereas 
the reversed order is calculated for MgB2 (see Table ||) 
although both compounds are structurally quite simi- 
lar. This effect was pointed out before and .studied in 
Mgi_2;Al x B2 mixed crystals by Renker et aZ.j£j The in- 
terchange, which occurs only in undoped or moderately 
substituted material (0 < x < 0.2), was also attributed 
to the electronic changes, especially the disappearance of 
the hole pockets from the Fermi surface for x > 0.2. 

The metal content in MgB2 affects the band filling, a 
larger Mg deficiency moving the Fermi level further be- 
low the top of the a bands. It is therefore to be expected 
that the anharmonicity of the Ei g mode should decrease 
with decreasing Mg content. Lattice dynamical calcu- 
lations showed that the— E-2 g anharmonicity decreases 
with increasing pressuretj and the initially very large 
mode-Griineisen parameter decreases too.Ea If the mode- 
Griineisen parameters decreases as function of band fill- 
ing at ambient pressure, i.e., due to non-stoichiometry of 
MgB 2 , it would, qualitatively, lead to the observed re- 
lation between dT c /dP and Mg deficiency. This effect 
would be a manifestation of the changes of the lattice 
dynamics rather than changes of the electronic density 
of states. A more detailed and quantitative analysis is 
certainly needed, but the present results are indication of 
the importance of stoichiometry with regard to the super- 
conducting properties of MgB2 , specifically the pressure 
dependence of T c . 



V. CONCLUSIONS 

We have studied the crystal structure of AIB2 by x- 
ray powder diffraction to 40 GPa. The compressibility 
is moderately anisotropic consistent with the anisotropic 
bonding properties. In the pressure range studied here 
we did not observe a structural phase transition. Our 
x-ray diffraction data indicate an Al deficiency of ^11% 
in agreement with previous reports. Despite the neglect 
of this non-stoichiometry in our first-principles calcula- 
tions, the calculated structural properties are in good 
agreement with the experiment. 

The Ei g zone-center phonon in metallic AIB2 can be 
observed as a well-defined Raman peak. We conclude 
that the lack of such a Raman feature in MgB2 is neither 
related to the metallicity or disorder of the sample nor 
is it a generic property of AlB2-type compounds. Our 
observations rather support the view that it is due to the 
strong electron-phonon coupling and/or anharmonicity 
which are distinct properties of MgB2- We found some 
deviation of the calculated pressure dependence of the 
Ei g phonon frequency of AIB2 from the experimental 
data and tentatively attributed this to the Al deficiency 
of the AIB2 sample which was not taken into account in 
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the theory. 

Correlations between non-stoichiometry of MgB 2 and 
its superconducting properties have been pointed out pre- 
viously. Here we considered possible effects of Mg defi- 
ciency in MgB2 on its electronic structure and lattice dy- 
namics. The anticipated changes are consistent with the 
available experimental data on the correlation between 
Mg content and the pressure dependence of T c . This 
leads us to propose that the large variation of the pres- 
sure dependence of T c (-0.7 to -2.0 K/GPa) in MgB 2 
in various experiments may be caused by the effect of 



non-stoichiometry on the lattice dynamics, mediated via 
changes in the electronic structure of MgB 2 . 
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